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Abstract 

This paper provides asymptotic structure at spatial infinity of plane steady 
Stokes flow in exterior domains when the obstacle is rotating with constant 
angular velocity. The result shows that there is no longer Stokes’ paradox due 
to the rotating effect. 


1 Introduction 

Let be an exterior domain in the plane M 2 with smooth boundary dfl, and consider 
the motion of a viscous incompressible fluid around an obstacle (rigid body) K 2 \ D. 
As compared with 3D problem, we have less knowledge about exterior steady flows in 
2D despite efforts of several authors mentioned below. The difficulty is to analyze 
the asymptotic behavior of the flow at infinity. This is related to the following 
hydro dynamical paradox found by Stokes (1851): The problem 


— A u + Vp = 

0, div u — 0 in O, 

(i.i) 

o' 

cs 

CO 

?3 

u —» u Q o as x —» oo 

(1.2) 


admits no solution unless u a0 = 0, where u(x) = (ui,U 2 ) t and p(x) denote the 
velocity and pressure, respectively, of the fluid. Throughout this paper, all vectors 
are column ones and (-) T denotes the transpose of vectors or matrices. Later on, 

* Dedicated to Professor Reinhard Farwig on his 60th birthday 


1 



Chang and Finn [6] made it clear that the Stokes paradox is interpreted in terms of 
the total net force exerted by the fluid to the obstacle 

N:= T(u,p)u da, (1.3) 

J dfl 

where T{u,p) is the Cauchy stress tensor given by 

T(u,p) \= (T jk (u,p)) = Du - pi, Du \=X7u+(X7 u) t , I=(5 jk ), (1.4) 

and v denotes the outward unit normal to dfl] in fact, they proved that the flow 
satisfying (11.111 can be bounded at infinity only if the net force (II ,3[) vanishes. This 
is an immediate consequence of asymptotic representation at infinity of solutions 
to (11.11) due to themselves [6], see (I2.12p . The original form of the Stokes paradox 
mentioned above follows from the result of Chang and Finn as a corollary because 
the net force (1 1.3 [) never vanishes provided that {u,p} is nontrivial and satisfies 
(11.11) together with u\qq = 0. There are some other forms of the Stokes paradox, 
see Galdi [US V.7], Kozono and Sohr [281, Theorem A]. 

For the case in which a constant velocity u Q 0 6 K 2 \ {0} is prescribed at infinity 
or, equivalently, the obstacle is translating with velocity —Moo (while the flow is 
at rest at infinity), Oseen (1910) proposed his linearization of the Navier-Stokes 
system around Moo to get around the Stokes paradox. This works well because the 
fundamental solution of the Oseen operator — Am+Moo-Vm+Vp possesses some decay 
structure with wake, while the Stokes fundamental solution grows logarithmically 
at infinity, see (12.131) . Finn and Smith [15] . [ 14] , [31] actually adopted the Oseen 
linearization to succeed in construction of the Navier-Stokes flow when m^ is not zero 
but small enough (and the external force is small, too, unless it is absent). Later 
on, Galdi na refined the result by means of ZA-estimates, see also m Section 
XII.5]. The similar existence theorem for the case m^ = 0 is still an open question 
even for small external force. Even before the results mentioned above, Leray [22] 
constructed at least one Navier-Stokes flow with finite Dirichlet integral without any 
smallness condition, however, the asymptotic behavior at infinity of his solution is 
still unclear and all the related results obtained so far are partial answers (Gilbarg 
and Weinberger |2T] , [22] and Amick [T], [2] ). For details, see Galdi [IS. [18| and [T9] . 
It should be noted that symmetry helps to attain the boundary condition u —> 0 at 
infinity, see m , m , m and the references therein. Among them, Yamazaki [22] 
employed a linearization method to construct a small Navier-Stokes flow decaying 
like l^l" 1 at infinity under a sort of symmetry; indeed, the symmetry he adopted 
enables us to avoid the Stokes paradox since the net force vanishes. 

In this paper it is shown that, instead of the translation mentioned above, the 
rotation of the obstacle leads to the resolution of the Stokes paradox in the sense 
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that: (i) The flow can be bounded (and even goes to a constant vector at the 
rate |x| _1 ) at infinity even if the net force (11.31) does not vanish (Theorem 12.11) : 
(ii) Given external force decaying sufficiently fast, there exists a linear flow which 
enjoys u{x) = 0(|:c| -1 ) as |x| —* oo (Theorem 12.21) . We also provide a remarkable 
asymptotic representation of the flow at infinity, in which the leading term involves 
the rotational profile x L /\x\ 2 whose coefficient is given by the torque, where x L = 
(— X 2 ,Xi) t , see (I2.10|) and (12.161) . Here, the linear system arising from the flow 
around a rotating obstacle with constant angular velocity a G R \ {0} is described 
as 

— A u — a (x -1 • Vm — n -1 ) + Vp = /, div u — 0 in Q (1.5) 

in the reference frame attached to the obstacle. We recall the derivation of (II. 5ft in 
the next section. 

The essential reason why there is no longer Stokes’ paradox is asymptotic struc¬ 
ture of the fundamental solution of the system (II.5|) in the whole plane M 2 . Roughly 
speaking, the oscillation caused by rotation yields better decay structure of the fun¬ 
damental solution, from which combined with some cut-off techniques we obtain 
the main results. It is worth while comparing with the result [9] by Farwig and 
the present author on the 3D Stokes flow around a rotating obstacle, in which the 
axis of rotation (e 3 -axis without loss, where e 3 = ( 0 , 0 , 1 ) T ) plays an important role; 
in fact, 63 • N controls the rate of decay. The result would suggest better decay 
studied in this paper since we do not have the axis of rotation in 2D, however, there 
are some difficulties compared with 3D case. Look at (j3.2j) below, which would be 
heuristically a fundamental solution, but this is by no means trivial because of lack 
of absolute convergence unlike 3D case. We thus employ the centering technique 
due to Guenther and Thomann [23], that is, we subtract the worst part, whose 
time-integral converges on account of oscillation, from the integrand of (13.2(1 such 
that the remaining part converges absolutely and can be treated in a usual way. 
This technique is also needed to justify some estimates of the fundamental solution, 
see Lemma 13.31 Asymptotic analysis of the fundamental solution to find the asymp¬ 
totic representation (13.341) is similar to the argumant adopted for 3D ([9]), in which, 
however, the external force / is assumed to have a compact support. In this paper 
we will derive further properties of the fundamental solution and the corresponding 
volume potential (I3.65p to deal with the external force decaying sufficiently fast for 
|x| —» oo, see ( 12 . 6 p and ( 12 .9p . 

This paper is a step toward analysis of the Navier-Stokes flow around a rotating 
obstacle in the plane. To proceed to the nonlinear case, it is reasonable to consider 
the external force / = div F with F(x) = 0{\x\~ 2 ) in view of the nonlinear structure 
M 1 Vm = div (u®u), see Remark 13.21 This will be discussed in a forthcoming paper. 
As for asymptotic structure of the Navier-Stokes flow around a rotating obstacle in 
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3D, the leading term at infinity was found first by and then the estimate of the 
remainder was refined by [8]. 

When the rotating obstacle is the two-dimensional disk and the external force is 
absent, the Navier-Stokes system subject to the no-slip boundary condition (12.191) 
admits an explicit solution (12.20)) in the original frame, see Galdi [19], p.302]. Re¬ 
cently, Hillairet and Wittwer [25J considered small perturbation from this solution 
with large angular velocity |a| to find the Navier-Stokes flow decaying like \y \~ 1 at 
infinity, whose leading profile is given by y L /\y\ 2 . See also Guillod and Wittwer [24, 
Section 4], who provided among others numerical simulations of the related issue. 

This paper is organized as follows. In the next section, after recalling the equa¬ 
tions in the reference frame, we present the main theorems. Section 3 is essentially 
the central part of this paper and we carry out a detailed analysis of several asymp¬ 
totic properties of the fundamental solution of (11.51) in the whole plane K 2 . Section 
4 is devoted to decay structure of the system (11.51) to prove Theorem 12.11 In the 
final section we show the existence of a unique linear flow which goes to zero as 
|x| —» oo to prove Theorem 12.21 


2 Results 

We begin with introducing notation. Set B p (x o) = {x G M 2 ; \x — xq\ < p}, where 
Xo G K 2 and p > 0. Given exterior domain D with smooth boundary dfl, we fix 
R > 1 such that M 2 \ C Br( 0). For p > R we set = O D B p ( 0). Let D be one 
of f2, M 2 and Q p , and let 1 < q < oo. We denote by L q (D ) the usual Lebesgue space 
with norm || ■ ||l9(_d). It is also convenient to introduce the weak-L 2 space L 2,00 (D ) 
(one of the Lorentz spaces, see [3]) by L 2,00 (D ) = (L 1 (D), L°°(D)) 1 j 2 ^ with norm 
II • \\l 2 ’°°(d), where (■, -)i/ 2 ,oo stands for the real interpolation functor. The measurable 
function / belongs to L 2, °°(D ) if and only if sup T>0 r|{x G D\ |/(x)| > t}\ 1 ^ 2 < oo, 
where | ■ | denotes the Lebesgue measure. Note that L 2 (D ) C L 2,0O (D); indeed, 
|x| _1 G L 2 ’°°(D). By H k (D) and Hq(D) we respectively denote the L 2 -Sobolev 
space of k -th order (k > 1) with norm || • \\u k {D) an d fl ie completion of C£°(D) 
(consisting of smooth functions with compact support) in H l (yD). We use the same 
symbol for denoting the spaces of scalar, vector and tensor valued functions. 

Before stating our results, we briefly explain the derivation of the system (11.51) 
for the readers’ convenience although that is the same as in 3D case (53, EES). 
Suppose a compact obstacle (rigid body) K 2 \ D is rotating about the origin in 
the plane with constant angular velocity a G M \ {0}, and let us start with the 
nonstationary Navier-Stokes system 

d t v + v ■ Wv = Av — Vg + g, div v — 0 
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in the time-dependent exterior domain hl(t) = {y = 0(at)x ; x G hi} with 


0(t) = 


cos t — sin t 
sin t cos t 


where v(y,t ) and q(y,t) are unknowns, while g(y,t ) is a given external force. The 
fluid velocity is assumed to attain the rotational velocity of the rigid body on the 
boundary <9hl(f) (no-slip condition), while it is at rest at infinity, that is, 

v|an(t) = ay L , v -> 0 as \y\ ->• oo. 

We take the frame attached to the obstacle by making change of variables 

y = 0(at)x , u(x, t) = 0(at) T v(y, t ), p(x, t) = q(y , f), 
f(x,t) = 0(at) T g(y,t), 


so that the equation of momentum is reduced to 

d t u = 0(at) T d t v + 0(at) T ^aO(at)xj • V,/c + aO(at) T v 

= 0(at) T (—v ■ X7 y v + A y v — X7 y q + g) + a (x -1 • X7 x u — u ± ) 

= —u ■ \7 x u + A x u — V a .p + f + a (x -1 • V x u — u L ) 

in hi, where 0(f) = ^O(f). If / is independent of f, then one can consider the steady 
problem 


— Au — a (or 1 • Vu — u ± ) + Vp + u ■ Vu = /, div u — 0 in hi (2.2) 


subject to 

u |an = ax -1 , u —> 0 as |x| —> oo. 

It is sometimes convenient to write the LHS of (12.21^ as divergence form 

A u + a (x 1- ■ Vu — w 1 ) — Vp — u ■ V« 

= div (Siu.p) - u® u) = I { S jk (u,p ) - UjU k } 


(2.3) 


vfc=l 


3 = 1-2 


with 

S(u,p) = (, S jk (u,p )) = T(u,p ) + a (u ® x 1 - — x -1 (8) u) ( 2 . 4 ) 

where T(u,p) is given by ( 11 . 41 ) , u®v — ( UjVk ) stands for the matrix for given vector 
fields u and u, and <9*, = d Xk . 
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The only problem we intend to address in this paper is the associated linear 
system (11.50 . On account of the relation 


/ [(x' L ■ \/u — w 1 ) ■ v + u ■ (x ± ■ Wv — v 1 )] dx = / [y ■ x _l )(m • v) da (2.5) 
Jn Jan 

for vector fields u and v (so long as the calculation (14. 12)1 in section 0] makes sense), 
the opeartor u ^ x L ■ Vu — w 1 is skew-symmetric under the homogeneous boundary 
condition. Also, by using the auxiliary function (j5. 1 j) below, our problem with 
boundary condition (12..311 L can be reduced to the problem with the homogeneous one. 
Hence, it is not hard to find at least one solution with Vw G L 2 {Vt) for (11.51) (even for 
the Navier-Stokes system (12.2)1 without restriction on the size of |a|) subject to the 
boundary condition u\on — ax ± (only (12.3p i ) along the same procedure as in Leray 
[2Sj provided / = div F with F G L 2 (fl), however, we do not know whether the 
behavior (12.31) 2 at infinity is verified. The asymptotic structure of this solution for 
/ decaying sufficiently fast at infinity and, more generally, that of {u,p} satisfying 
(11.51) without assuming any boundary condition on dfl are given by the following 
theorem. For simplicity we are concerned with smooth solutions although the result 
can be extended to less regular solutions (in view of Proposition 13.21 for the whole 
plane problem). 

Theorem 2.1. Let a G M\{0}. Suppose that {u,p} G Hl oc (fl) x L 2 oc [Ll) is a smooth 
solution to the system (11.51) with f G C°°(r2) satisfying 


\x\\f(x)\dx < oo, |/(x)| < 


C 


(1 + |x| 3 )(log (e + |a;|)) 


where the constant C > 0 is independent of x G Q. Assume either 

(i) Vw G L r (Ll \ B r ( 0)) for some r G (1, oo) 
or 

(ii) ufx) = o(|x|) as |x| —* oo. 

Then there are constants u ^ G M 2 and p ^ G M such that: 


( 2 . 6 ) 


1. (asymptotic behavior) 


u(x) = Uoo + (1 + |a| x ) 0(|a:| x ), 
p(x) = -au^ ■ x + Poo + 0(|x| _1 ), 


as |x| —>■ oo. 

2. (energy balance) 
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We have Va 6 L 2 (Q) (even if we do not assume (i) with r — 2) and 


- / \Du\ 2 dx = 


n \ / s aiv-x 1 -). 

( T{u,p)u) ■ ('U - Moo) H --- 1 U - Moo 

ion L 2, 

+ / / • (m - Moo) dx 


da 


( 2 . 8 ) 


with T(u,p ) := T(u, p + au ^ • x — Poo), where Du and T(-, •) are as in (11.41) . 


3. (asymptotic representation) 

If in addition 

f(x) = o(|x|^ 3 (log |x|) _1 ) as |x| —> oo, 

then 


u{x ) Mqo 


ax 1 — 2 fix 
Att\x \ 2 


+ (1 + |a| 3 ) o(|x| x ) 


as |x| —> oo, 


where 


a= y 1 ■ {(T(u,p) + om® y 2 -) z/} da y + / y ± ■ f dy , 
Jdo Jo 


/3 — v ■ m da. 

JdVL 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


If in particular f G C^°(f2) ; that is, the support of f is compact in M 2 , then 
the remainder decays like 0 (|:r| _2 ) m ( 12 . 101 ) . 


Note that T{u,p)v belongs to H~ l / 2 (dD) := H l ^ 2 {dD)* by the normal trace 
theorem since T(u,p) G L 2 (Qr) and div T(u,p ) = —a (or 1 • Vn — m - 1 ) — / G 
L 2 (h2/{). Therefore, the boundary integral f dn y L • (T(u,p)u ) da y can be understood 
as T(u,p)v) H -i/ 2 ( d n' ) in (12.111) . Since u G i7 1//2 (<9Q), the same reasoning 

as above justifies f dn ( T(u , p)z/) • (m — Mqo) da in (12.81) . All the other integrals in (12.81) 
and (12.111) also make sense. 

For the usual Stokes system (11.11) . under the same growth condition on u(x) as in 
Theorem 12.11 there is a constant m^ G M 2 such that (Chang and Finn [6j Theorem 
1 ]) 

u(x) = Moo + E(x)N + 0(\x\ 3 ) (2.12) 

as |x| —> oo, where 

£(*) = fz 


, 1 \ , x <g) x 

log —7 I + 


\x\ 


\x\ 


(2.13) 
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is the Stokes fundamental solution and N denotes the net force (II.dp . We observe 
the remarkable difference between (12.71) , and (12. 12j) : in fact, the flow is bounded in 
Theorem 12.11 even though the net force N does not vanish. We would say that this 
is the resolution of the Stokes paradox. 

The leading term of (12.101) is of interest since it contains the rotational profile 
x^/\x\ 2 , which comes from the leading term of the fundamental solution of (11.51) . 
see (I3.34p . The other profile — x /{ 2 ' k \ x \ 2 ) is called the flux carrier. If in particular 
the flux /3 at the boundary vanishes, then the leading term is purely rotational and 
that is just the case in the next theorem. Look at the coefficient (12. lip of x^/\x\ 2 , 
where the integral fg^y 1 ■ (T(u,p)v ) da y stands for the torque exerted by the fluid 
to the obstacle. It is worth while noting that, in three dimensions, one finds the 
rotational profile (e 3 x x)/\x\ z , whose coefficient involves the torque, in the second 
term after the leading one. For details, see Farwig and Hishida P Theorem 1.1], It 
is reasonable that both x 1 -/\x\ 2 = V x log |x| and x/\x\ 2 = V log |x| are solutions to 
(11.51) with / = 0 in M 2 \ {0} together with the constant pressure and, therefore, so 
is the leading term of (12.101) . In fact, we observe 

A—= 0, ^ = V4-, div nj = 0 m K 2 \ {0} 

as well as (14.ip (with x 0 = 0) below. 

In Theorem 12.11 it is also possible to find the asymptotic representation of the 
pressure p(x) without assuming any growth condition on p(x) itself since it can be 
controlled by the growth of u(x) via the equation (11. 5p ,. The leading profile of 
p(x) + aw x ■ x — in (I2.7P 2 is just the fundamental solution Q(x) = 2 ^\ x \ 2 the 
pressure to the Stokes system. This is because 

div (x x • Vw — u x ) = x 1 - ■ Vdiv u = 0 (2.14) 

so that the pressure part of the fundamental solution is independent of a € M. Thus 
we are not interested in the asymptotic representation of the pressure, which the 
rotation of the obstacle does not affect so much. The coefficient of the leading profile 
Q(x) is rather complicated in Theorem 12 . II but it becomes just the force in the next 
theorem, see (12.171) . 

The next question is whether one can actually construct a solution to (| 1.5 p when 
zero velocity is prescribed at infinity as in (12.3p . The following theorem gives an 
affirmative answer. 

Theorem 2.2. Let a G M\{0}. Suppose that f = div F E C°°(Ll), with F E L 2 {Vt), 
satisfies (12.6p . Then the system (II. 5p subject to (12. 3 p admits a smooth solution 
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{u,p}, which is of class u G L 2, °°(f2) D H( nc (Q), p G L 2 oc (Vl) as well as Vu G L 2 (H) 
and fulfills 

1 + M + (1 + M *) (ll-^1U 2 (fi) 

+ Jj x \\f( x )\dx + sup \ x \ 3 (fog (e+\ x \))\f( x )f) , ( 215 ) 

IIVm||l2(q) < ||F|| L 2 (n) + C|a|, 

with some C > 0 independent of f and a G M \ {0} ; and {u,p} enjoys all the 
assertions in Theorem \2.1\ with {woo,Poo} = {0,0}. In particular, we have 

u{x)= ( [ y 1 - ■ (T(u,p)v)d<r y + [ y^ ■ f dy\ X |? + (l + |a|~ 1 )o(|a;|~ 1 ), (2.16) 
\Jdil JQ J 47r Fl 

P( x )=([ T(u,p)uda+ [ fdy] • X +0(|x|~ 2 ), (2.17) 

\Jdci Jn J *n\x\ 

as |x| —> oo under the additional condition (12.91) (which is needed only for (12.161) I. 
This is the only solution in the class Vu G L 2 (Q), {u,p} G L 2 oc (fiX) with {u,p} —$■ 
{0, 0} as |x| —> oo. 

Note that, when a = 0, the problem is not always solvable for given external force 
/ = div F even if F G C{fi (hi), that may be also regarded as the Stokes paradox. 
The L^-norm of |x| \u(x) | away from the boundary can be also estimated by the RHS 
of (12.15j) x (see (15.61) for an approximate solution). In order to control the L^-norm 
of u(x) near the boundary dQ, the class H( oc (Q) is not enough. We put the term 
x 1 - ■ Vu — u 1 - in the RHS and use the regularity theory of the usual Stokes system 
up to the boundary to show that u G H( oc (Q) C f!) together with a certain 
estimate, which enables us to obtain the similar estimate of sup xGn (l + |x|)|w(x)| to 
(l2T5D r 

We conclude this section with the following exact solutions of both the Stokes 
and Navier-Stokes boundary value problems without external force. The Stokes 
flow (I2.21|) seems to be well known since it is found in some old literature. The 
Navier-Stokes flow (I2.20|) is found in the second edition of [19], p.302] (I learned it 
from Professor Masao Yamazaki around 2008). Suppose the unit disk (rigid body) 
Ri(0) is rotating about the origin with constant angular velocity a G K \ {0}. Then 
the Navier-Stokes flow in the exterior H = M 2 \ -£^(0) obeys 

— An + Vg + v ■ Vu = 0, div v — 0 in hi (2-18) 


|w||l2,c 


yo) < C 
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subject to 

f|ao = ay- 1 , v -P 0 as \y\ ->■ oo 
and this problem has a solution 

I 2 

/ \ oy / x ~ a 

v (V) = QW = + constant. 

Also, the associated Stokes problem 

—An + Vg = 0, div v — 0 in 12, 
sibject to (12.191) admits a solution 

ay 1 - 


v(y) = 


\y \ 2 ’ 


q(y) = constant. 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


Note that the Stokes flow (I2.2ip does not contradict the Stokes paradox because 
f dQ T(v, q)v do = 0 due to symmetry. When the obstacle is a disk, we do not 
necessarily have to make the change of variables (12.ip . nevertheless, we can do so 
and this case is not excluded in the present paper. Steady flows in the original frame 
correspond to time-periodic flows and are not steady in general in the reference frame 
via (E3J. But the Stokes flow (12.211) becomes the steady one u{x) = ax 1 -/ |x| 2 , p(x) = 
constant in the reference frame as well and this may be regarded as a special case 
in Theorems 12.11 and 12.21 (when we take p — 0 in the latter theorem); indeed, one 
can verify 

/ y _L • ( T(u,p)v ) d<T y = 47ra 
Jon 

in (12.161) . Recently, Hillairet and Wittwer [25j proved that if the boundary value 
u|a £2 is sufficiently close to ay 1 - with |a| > \/48 in a sense and f dQ v ■ v da = 0, then 
the Navier-Stokes system (12.181) in the exterior 12 = M 2 \ -Bi(O) of the disk subject 
to this boundary condition admits at least one smooth solution, which decays like 
\y\~ l as \y\ —» oo. The leading profile of their solution is given by y _L /||/| 2 with some 
coefficient close to a. 


3 Fundamental solution 


In this section we derive the decay structure of the fundamental solution of the 
linear system (1 1.5 p in the whole plane M 2 when a G R \ {0}. Because of (12.141) the 
pressure part of the fundamental solution is 


Q(x - y) 


x-y 

2n\x — y\ 2) 


(3.1) 
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while the velocity part is given by 


r a (x,y) 


0(at) T K(0(at)x — y, t) dt, 


(3.2) 


where 

K(x, t ) = G(x, £)I + H(x , i) 

is the fundamental solution of unsteady Stokes system (a = 0), and it consists of 
the 2D heat kernel 

G(M) = _L e -MV« 

and 2x2 matrix 


H(x, t ) 


V 2 G(:r, s) ds 



x ® x 
4s 2 



(3.3) 


In 2D case one can write (13.3jl in terms of elementary functions 


H(x, t ) 


=tlp. G{xJ) + 

\xr 


f x® x I \ 1 — e hi 2 / 4 * 

\ \x\ 2 2) 7i\x\ 2 


(3.4) 


while one cannot in 3D, see [9|. One needs more careful argument than 3D case p] 
to prove that (13.2p is actually the fundamental solution, see Proposition 13.21 

Indeed the integral representation (13.211 does not absolutely converge, but it is 
convergent due to oscillation 0(at) T with a € R \ {0}, see Lemma 13.21 This is 
a contrast to the case a = 0, in which (13.211 is not convergent. In this case one 
needs the centering technique to recover the convergence, which leads to the Stokes 
fundamental solution E(x) given by (I2.13P as follows: 



K(x, t ) 


P - e/ 4 * \ 

--I) 

8nt J 


dt = E(x). 


(3.5) 


This was clarified by Guenther and Thomann [23) Proposition 2.2], As a part of this 
technique (13.511 . the fundamental solution of the Laplace operator in two dimensions 
is recovered exactly as 





(3.6) 


in terms of the heat kernel, see m Lemma 2.1]. Although we do not need the 
centering technique in the representation (13.211 itself, we will use this technique to 
justify some formulae in this section. 


11 


















Remark 3.1. In [TT, . p.301] Farwig, Hishida and Muller mentioned that the integral 
kernel J 0 °° 0{t) T G(0(t)x — y,t ) dt should be modified to recover the convergence in 
two dimensions. But this is redundant as we will see in Lemma \-i. 21 by making use 
of the oscillation. 

For convenience we will collect a few elementary fonmnlae, which will be used 
several times. We omit the proof that is nothing but integration by parts. In the 
first assertion below it is possible to derive even faster decay r -{ 2 ( m ^ 1 )+ 2fc l for every 
k G N by fc-times integration by parts, but (13.71) and (13.81) are enough for later use. 
Note that they are not absolutely convergent for m < 1 (the only case we need is 
m — 1 ). 

Lemma 3.1. Let r > 0. 


1. Let a G M \ {0} and m > 0. Then 


e iab e -r*/t &1 


< 


C 


\a\r 


2 m ' 


coo coo 

e iat / e -r 2 /s 


. 2 /. ds 


dt 


< 


C 


2m - 


/ 0 J t s m+1 \a\r 

with some C = C(m ) > 0 independent of r > 0 and a G 

i = x/—T- 


(3.7) 

(3.8) 
\ {0}, where 


2. Let m > 1. Then 



7 (m — 1 ) 

rp2(m— 1 ) ’ 



7 (m — 1 ) 

,j.2(rn — 1 ) ’ 


where 7 (-) denotes the Euler gamma function. 


(3.9) 

(3.10) 


We begin with the following lemma, from which the function (13.2|) is well-defined. 

Lemma 3.2. Let a G M \ {0}. Then the integral T a fx,y) given by (13. 2 p converges 
for every ( x , y) Gl 2 x M 2 with x j- y. 
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Proof. We decompose T a (x,y) as 

r„(z, y) = rsti, y) + p(x, y ), ri(i, y ) = r“(i, i/) + P 2 (*,</) 

with 


r°(x, y)= / 0(at) T G(0(at)x — y,t) dt, 


r a( x ^y)= / °( at ) / G(0(at)a;-y,s) 


( 0{at)x — y) ® ( 0(at)x — y) 


4s 2 


ds dt , 


poo /»oo 1 

\T 1 ' 1 


( 2 ;,?/)= / O(at) 1 / G(0(at)x — y, s)——ds dt. 

Jn Ji 


(3.11) 

We start with the convergence of T°(x, ?/) by using the centering technique as in 
(13.(jp . By (13.7p we know that 


O(at) e 


T-l/At 


dt 


< 


C 


Hence, it suffices to show the convergence of 

J 0(at) T (^ e -\ 0 (at)x-y\pU _ e -l/4^ ^ 

As we will see, this is absolutely convergent. For large t, we have 


(3.12) 


(3.13) 


D -\0(at)x-y\ 2 /4:t _ -1/4 t 


dt 


dt 


J < \\0{at)x-y\ 2 -l\ — 


< 


(\x\ + \y\Y + l 


For small t, we use the relation 

1 0(at)x — y | 2 — \x — y | 2 + 2 at[0(a.9t)x) ■ ( 0{a6t)x — y) 
for some 9 = 6(a,t,x,y ) G (0,1), where Off) = ^ 0(t ). Then we have 


e —|0(afp—?/| 2 /4t _ g—1/41 


dt 


< 



Hz-2/| 2 /4t e MM(M+M)/ 2 _|_ g -l/4t 


dt 


(3.14) 


(3.15) 


with 


e -|x-y| 2 /4t ^ = 

t 


1 0 Jl 


<4 + 


'■Vh-d 2 dt 
t 


= 4 + 2 log 


•i/h-?/l 2 ^ 

| —l/4( ^ 

) 

1 t 

1 

(0 < £ 

x - 2/| 

13 



(3.16) 













while 

[ e -\x-y?/*t & < 4 (|a? — 3/| > 1). 

JO t 

This concludes that (j3.13|) is absolutely convergent. 

The next integral T A1 (x,y) is absolutely convergent without centering as above. 
Given (x, y ) with x j- y, there is 5 = 5 (a, x, y) > 0 such that 

0< |3; ~ y| <\0(at)x-y\ 2 < ^ X ~ y \ , 0 < Vt < 5, (3.17) 

on account of lim 1 0(at)x — y I 2 — \x — y I 2 . This together with (13.1011 implies that 
t->o 


poo poo 

L l 1 

rS roo 


r \0(at) X -y\z/4s \Q(at)x - V \ 2 ^ ^ 


< 


< 



0 Jt 


e -\x- y \ 2 /&s 3 l x y \ 2 dsdt + 

2s 3 


poo poo 


'5 Jt 


(M + M) 5 


ds dt 


3k-»i 2 r r _- b -,iv8, ds dt| (m + m) 2 rat 


(3.18) 


o Jt 
2 


t 2 


= C + 


(M + \y\) 

25 


Finally, similarly to the argument of convergence of T®(x,y), we can discuss 


T 12 {x,y). From (13.811 it follows that 


ds 


C 

< — 

a 


0(at) T / e l / As — dt 

o Jt s 

It thus remains to show the convergence of 

r 0{at) T r (e-\0{at)x-y\>/As _ g-1/4^ J dt 


(3.19) 


(3.20) 


For large t, we have 


POO poo 


n Jt 


e -\0(at)x-y\ 2 /4s _ g -l/4s 


ds , 
-jdt < 
S Z 


< 


ds 


\\0(at)x — y\ — 1 —— dt 


(N+W + i 


4s 3 
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For small t, as in f)3.15j) . we use (13. 14j) to find 


with 


1 [‘ OO 



' 0 Jt 

noo 


e -\0(at)x-y\ 2 /4s _ g-1/4 s 


ds 


dt 


< 



0 Jt 


e -\x-y\ 2 /is g|o||x|t(|s|+|j/|)/2s _j_ g-l/4s 


ds 


dt 


< e \ a \\x\{\x\+\y\)/2 


1 /*oo 



0 Jt 


e -i,-„p/4» if + 4 


a OO 


e -\x-y\ 2 /is < l±dt= ( / + 

ol/\x-y\ 2 


1 rl/\x—y\ 2 ' 

+ / 

0 J 1 


e -i/ 4 s dt 
s z 


< 4 + 4 / (1 - e -1/4s ) ds < 4 + 


s 


(3.21) 


(3.22) 


= 4 + 2 log 


while 


x-y | 


1 /»oo 



(0 < |x - y\ < 1), 


e -\x-y\ 2 /4s ds dt < 4 (|x — ?/!>!). 


0 Jt 


This implies the absolute convergence of (13.20)1 . We have completed the proof. □ 


We have concentrated ourselves only on the convergence of the integral (13.2)1 . So 
the estimates appeared in the proof above are not related to the asymptotic behavior 
with respect to (x, y ) at large distance, which will be discussed in a different way in 
Proposition 13. 11 but we have tried to derive the singular behavior for \x — y\ —y 0 as 
less as possible, see (13.16)) . (13.18)) and (13.22)) . This behavior should be logarithmic, 
otherwise (13. 2 1) cannot be the fundamental solution, but the behavior (13.18)) is not 
clear since 5 depends on x,y (probably, the part T 41 (a;,|/) would be bounded for 
\x — y\ —y 0 as in the second term of the Stokes fundamental solution (12.131) ). In 
order to ensure that the volume potential (13.65P below is well-defined, we will show 
the following lemma. The growth rate (13.23)) with p — 2\x\ will be also used to show 
asymptotic representation (13. 70)) below. 


Lemma 3.3. Let a G R\{0} There is a constane C > 0 independent of a G M\{0} 
such that 


\T a (x,y)\dy < C\a 


-1 „2 


P +Cp log p, 


\y\<p 


(3.23) 



W x T a (x,y)\dy < Cp, 


(3.24) 
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for every igK 2 and p > |x| + e. 


Proof. To this end, it is convenient to use another representation (I3.4[) of H(x,t ) 
together with the centering technique (13.5ft due to Guenther and Thomann [23]. But 
we subtract /8tt t) I instead of (e _e / 4t /87rt) I since there is no need to derive 

E(x). First of all, it follows from (13.71) that 


r°° P -i/4t 

/ 0(at) T - dt 

/o 1 ' 8vr t 


C 

< —. 
a 


(3.25) 


We set 


p -l/At 

r a{x,y) := I 0(aty \^K(0(at)x-y,t) - 


I dt. 


(3.26) 


We will see that the integral of this integrand over (0, oo) x B p { 0) with respect to 
(t,y) is absolutely convergent. By the transformation y = 0(at)z we have 


roc 

l y 

y\<p 

e" 1/4i 

K(0(at)x y,t ) g ^ I 

/»oo 

Jo J 

( 

\ z \<p 

p-l/4 1 

K (x — z,t) - 1 

K ’ 8-Kt 

dz dt 


dy dt 


(3.27) 


The useful decomposition discovered by [23j is 


-l/4t 

e -|x| 2 /4 1 _ g-1/4 1 / g-|x| 2 /4 1 ^ _ e ~\x\ 2 /&t 


8nt 


-1 + 


Then we find 


8nt 2ir\x\ 


\A\ dt = C I log |x| 


2x ® x 


\x\ 


A + B 


8n 


(3.28) 

(3.29) 


while we see from the transformation r = |x| 2 /4t that 


I B I dt — 


2x ® x 


\x\ 


—re T + 1 — e 


r z 


dr = 


2x ® x 


\x\ 


< C (3.30) 


since 


0 < 


-re T + 1 — e 


d f e 
dr 
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for every r > 0. By the Fubini theorem we obtain 



T a (x,y)\dy <C 



(1+ | log I# — 7/| |) dy 


<Cp 2 + C / | log \x — y \| dy 

J\y-x\<p+\x\ 

< Cp 2 + Cp 2 log p 


for p> \x\ + e. This together with (13.25(1 concludes (13.231) . 

For the estimate of V x r a (x, y), we first need to justify 

POO 

VjT a (x,y)= / 0(at) T S7 x [K(0(at)x - y,t)] dt (3.31) 

Jo 


with the aid of r a (x, y) given by (13.261) . Given p e Cg°(K 2 ) arbitrarily, we have 


(r o (-, y), div ip) 



K(0(at)x 


y,t ) 


8nt J 


div p 


dt 


because this integral over (0, oo) x B L ( 0) with respect to (t,x) is absolutely conver¬ 
gent by the same reasoning as in the proof of (13.23(1 . where L > 0 is taken in such 
a way that Supp p C Bl( 0). We then use 

/ OO 

s -5/2 e -p| 2 /ie S ds (3.32) 

together with (j3.9p - (j3.10p to get the absolute convergence 




W x [K(0(at)x — y,t)] | dxdt < C 




\(WK)(0(at)x 


y,t) | dx dt 


= C 

< c 


POO P 

/ / I 

Jo J\x\<L 



dx 

x-y | 


y,t )| dxdt 


as in (13.271) . Hence we obtain 

POO 

(r o (-,Z/),div p) = - / <0(at) T V a; [/l(0(at)x - y,t)],<p) dt 

POO 

0(at) T 'V x [K(0(at)x — y,t)] dt,p 


17 
















for all p G C^°(R 2 ), which implies (I3.3ip since W x T a (x,y) = W x T a (x,y). Once we 
have that, by the same reasoning as above we get 


[ \\7 x T a (x,y)\dy < C f f \(S7 K)(0(at)x - y,t)\dt dy < C ( dV 

'\y\<p J\y\<pJo J\y\<p\ x d I 


which leads to (13.24[) for p > |x| + e. 


□ 


The following estimate provides the decay structure of T a (x, y) and plays a crucial 
role in this paper. 

Proposition 3.1. Let a G R \ {0}. 

1. There is a constant C > 0 independent of a G R \ {0} such that 

C'(|a|- 1 +M 2 ) 


r a(x,y) - 


x <8) y 


47t|x| 


< 


\x\ 


(3.33) 


for all (x,y) G K 2 x l 2 with |x| > 2\y\. In particular, we have 


r a (x,y) = 


x ® y ,_ 2 


47t|x| 


+o{ i*r), 


(3.34) 


as |x| —y oo so long as \y\ < p, where p > 0 is fixed. 

2. Similarly, there is a constant C > 0 independent of a G R \ {0} such that 

Cda ]^ 1 + |x| 2 ) 


r a(x,y) 


x x (8) y 1 - 


47t| 2/p 


< 


M s 


(3.35) 


for all ( x , y) G R 2 x R 2 with \y\ > 2\x\ 


Proof. Tha latter assertion follows from the former one because T a (x, y) = T_ a (?/, x) T 
and {y L QPx- 1 ) 1 = x L ®y x . We will show (I3.33p . which immediately implies (I3.34p . 
Let us start with T°(x, y) given by (13.lip . We use the Taylor formula with respect 
to y around y = 0 to see that there is 6 = 9(a, t, x, y) G (0,1) satisfying 


e -\0(at)x-y\ 2 /4t 

_ r -\x\ 2 /4t , -\x\ 2 /4t ( 0(at)x ) • y 

+ 2 1 

, 1 —\o(at)x— 9 y\ 2 / 4 t t ( 0(at)x - Qy) ® (0(at)x - 6y) - 2 1 1 
~2 6 V At 2 y • 


(3.36) 
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According to this formula, we decompose V®(x,y) as 


r°(x,y) = r° 1 (o;) + r° 2 (x, 1/ ) + rf( a ;, 1 /). 

It follows from (13.711 that 


1^)1 = 


1 /'°° 

4vr J 0 


0(at) T e~ lxl2/4t 


dt 


< 


C 


\a\\x\ 


Since 

and, thereby, 


(O(at)x) ■ y — {x ■ y) cos at + (or 1 ■ y) sin at 


(3.37) 

(3.38) 


{(O(at)a) ■y}0(at) T = - 


x ■ y 

x 1 - ■ y N 

) + 

cos 2 at j 

'x-y 

—x 1 - ■ y 

x-y / 

2 \ 

K X ± -y 

sin 2 at ( 

x 1 - ■ y 

x ■ 

■ y \ 


2 V 

-x ■ y 

X 1 - 

■y J ’ 



(3.39) 


we have 


C (£,!/) = 


,-|*l a /4t 


2 . i i dt 


167T j 
1 

47r|a:| 


x ■ y x • y 


t 2 \ —x ■ y x ■ y 

+ M° a 2 (x,y) 


+ *C(x,y) 


x ■ y x ■ y 
-x 1 - ■ y x ■ y 


(3.40) 


with 


\M° a 2 (x,y)\< 


C\y\ 


< 


C 


alld 3 laim 


(3.41) 


for |x| > 2\y\, which follows from (13.711 . Since e e v I 2 / 4 * < e M 2 / 164 f or |^| < 

\x\/2, it is easily seen that 


|r° 3 (x,r/)| <C\y\ 2 / (| x\H- 3 + t- 2 )e-W m dt 


C\yf 


\x\ 


(3.42) 


without using oscillation. Then (13.3711 . (13.4011 . (13.41 p and (j3.42|) imply that 


r2(*,y) 


x ■ y x • y 


47r\x\ 2 \ -x 1 - ■ y x ■ y 


<C{\a\- l + \y\ 2 )\x'- 2 


(3.43) 


for |x| > 21j/1. 
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We proceed to the decay structure of r„(x, y) given by (13. lip . Similarly to (13.36jl . 
we have the formula 



\0(at)x—y\ 2 /4s 

-\x\ 2 /4s i e -\x\ 2 /4s (O(at)x) ■ y 

2 S 

_|_ I e -\O(at)x-0y\ 2 /4s y T ( 0{at)x 


6y) ® (O(at)x 
As 2 


Qy) — 2s I 

- y 


(3.44) 


with some 6 = 9(a,t, s,x,y) E (0,1) and, correspondingly, we decompose r^x, ?/) 
given by (13.111) as 



x,y) = T a (hl/) + r/(i,t/) + r a (x,y). 


We write 

0(at) T [(0(at)x — y) ® ( 0(at)x — y)} 

— (x — 0(at) T y) <8) (O(at)x — y) 

. . , . . . . cos 2 at ~ sin 2 at ~ 

= A) + (cos at) A c + (sill at) A s H--— A c H--— A s 

with 


A 0 = Ao(x,y) 
A c = A c (x,y ) 

= A s (x, y) 
A c = A c (x,y ) 
As = A s (x,y) 


—3(x <g) ?/) + (a;- 1 <g> y- 1 ) 

2 ’ 

/ + J/l ®i®2 + 2/12/2 \ 

V X 1' X 2 + 2/12/2 + 2/1 / 


-X1X2 + 2/12/2 


+ 2/2 

^1^2 - 2/12/2 


-x • y 

x± ■ y ^ 

—x 1 - ■ y 

— X • V , 

—x -1 • y 

-x • y 

x • y 

—x -1 • y 




Using (13.8p and (13. 10p . we get 


r nl fa; y) = dll 

0 1 X,V) 16tt 


OO POO 



ds 


e -M As dt + M ni^ y) 


-3 {x®y) + (x ± ®y ± ) m 


87t|x| 


+ M a (U2/) 


(3.45) 


(3.46) 
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with 


(3.47) 


I M a ( x ,y)\< 


C 


ax 


for |x| > 2\y\. Look at ((31 


{( 0(at)x ) • y}0(at) T [(0(at)x — 2 /) <S) ( 0(at)x — y)\ = B 0 + (remainder) 


with 


„ x ■ y . x 1 - ■ y . x ■ y . , x L ■ y , ,. „ 

B 0 = H- 2 As = ~2^ X ® ^ H - 2 X 0 X + Bl = 


|x| 2 (x <g) y) 


+ B l 


which is independent of t, where B 1 is of degree one (resp. three) with respect to 
x (resp. y) and the remainder consists of all terms involving cos kat and sin kat 
(1 < k < 3). We thus find from (13.81) and (13.101) that 


r H*,j/) = 


|x| 2 (x <g y) 


OO poo 


647r 
x ®y 



e -w ■/‘•Li + C’hsl 


'0 Jt 

r!12 / 


47t|x 


2 + M a ( X , V ) 


with 


\M^(x,y)\< 


CiWl^ly] + bl 3 ) / C^a ]- 1 + \y\ 2 ) 


< 


\x\ 


\x\ 


for |x| > 2\y\. Without using oscillation, we see that 


00 poo 


K™(x,v)\<C\y\ 2 \x ' 2 



0 Jt 


(| X |V 5 + W 4 ) e-W 2/16s dsdt = 


C\y? 


\x\ 


for |x| > 2\y\. We collect (13.461) . (I3.47p . (I3.48p . (I3.49P and (I3.50p to hnd 

(x <g) y) + (x -1 <g) y A 


K^y)- 


87t|x| 


2 m „|-2 


< C'dal- 1 + | 2 /|^)|x 


for |x| > 2\y\. 

Finally, we decompose T^ 2 (x,y) given by (13.111) as 

r l 2 ( x , y) = r I 21 (x) + r 422 (x, y) + r^ 23 (x, y) 

by use of (I3.44p and deduce its decay structure. By (13.81) we have 

C 


|rr(*)l < 


ax 


2 ' 


(3.48) 


(3.49) 


(3.50) 


(3.51) 


(3.52) 
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As in the argument for T° 2 (o;, y), we employ (13.39P to obtain 


K 2 \x,y) = 


-1 

32^ 


r*oo poo 


e ~\x\ 2 /As ds dt l x-y X- ■ y \ + 122(x } 

' -x ■ y x-y 1 


f o Jt 


„_l_ 


87t|x| 


2 \ — 


with 


x-y x-y 
x 1 - ■ y x-y 


\Ml 22 (x,y)\ < 


(3.53) 


+ M 1 a 2 \x,y ) 


C\y I 


< 


C 


lallxl 3 |a||a ;| 2 


for |x| > 2\y\. Similarly to the argument for T 413 ^, y), it is seen that 

C\y \ 2 


r» oo poo 


T\ 2 \x,y)\<C\yY 


'0 Jt 


(|x|V 4 + s“ 3 ) e~W/ Ws dsdt = 


\x\ 


(3.54) 


(3.55) 


for |x| > 2\y\. We collect (13.5211 . (13.5311 . (13.54j) and (13.55)1 to obtain 

-1 ' — 


r ?(x,y)~ 


x■y x-y 
■x 1 - ■ y x-y 


87t|x| 2 

for |x| > 2\y\. Using the simple relation 

x-y x 1 - ■ y 


<C'(|a |- 1 + M 2 )|x | - 2 


(3.56) 


= x <S> y + x , 

-x^ ■ y x ■ y ' 

we gather (13.431) . (13.51 j) and (13.561) to conclude (13.34)) . The proof is complete. □ 

We next verify that ()3.2j) can be actually the fundamental solution. To this end, 
we need two lemmas. 

Lemma 3.4. Let f E L^R 2 ) n L°°(R 2 ) and 


p(x) = / Q(x-y) ■ .fly] dy, 


(3.57) 


where Q{x) is given by (13.11) . Set 


v°(x, t ) = 0(at) T / G(0(at)x - y, t)f(y ) dy , 


(3.58) 


v 1 (x, t ) = 0(at) T / H(0(at)x - y, t)f(y) dy , 


(3.59) 
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where H(x,t ) is given by (13.31) . Then they respectively satisfy 

d t v° + L a v° — 0, v°(;0) = f, (3.60) 

d t v 1 + L a v 1 = 0, u 1 (-,0) = -Vp, (3.61) 

in M 2 x (0, oo), where 

L a v := —Av — a ( x ± ■ Vv — . (3.62) 


Proof. The well-known estimate of singular integrals yields Vp G L q (J8L 2 ) for every 
q G (1, oo). By the derivation (12. ip of the equation (j2.2jl . it is obvious that v°(x, t ) is 
a solution to the Cauchy problem (13.6011 . where the initial condition is understood as 
lirn^o \\v°{t) — /||l<z(r 2 ) = 0 for every q G (1, oo). By the same reasoning, n 1 (x, t) = 
{v\,v\) T with 


v](x, t) = -'Y] 0(at) kj / G(0(at)x - y, t)d k p{y) dy 


= -/ d x .G(0(at)x-y,t )/ Q(y - z) • f{z) dz dy (j = 1,2) 


solves (13.6111 . Note that the integration by parts above can be justified since p G 
L r (M?) for every r G (2, oo) by the Hardy-Littlewood-Sobolev inequality. So we 
have only to deduce the representation (13.5911 . Using the relation 


Q{y) = 


y 


27t| y\‘ 


VG(y,r) dr 


and the semigroup property of the heat kernel, we find 

P P OO P 

ij(x,t)= V' / / d Xj G(0(at)x — y, t)(d m G)(y — z,t) dy dr f m (z) dz 

Jr 2 m Jo Jr 2 

P POO 

= - / y / d Xj d Zm G(0(at)x - z,t + r) dr f rn (z) dz 
Jm 2 m Jo 

P _ POO 

= iV' 0(at)kj / (d k d m G)(0(at)x - z,s) ds f m (z) dz 
J R 2 k,m J > 

which leads us to (13.5911 . 


□ 
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Lemma 3.5. Let e > 0. Let U G ^'(M 2 ) fulfill 

eU - AU -ax ± - S7U = 0 

where S' is the class of tempered distributions, 
denotes the Fourier transform of U. Similarly, if 

eu — A u — a(x _L • Vw — M" 1 ) + S/p = 0, 


Then Supp U C {0}, where U 
u G 5'(R 2 ) and p G 5'(K 2 ) satisfy 

div u — 0 in M 2 , (3.63) 


then Supp u C {0} and Supp p C {0}. 

Proof. We will prove the second assertion along the same idea as in EU, P3 Lemma 
4.2] (in which the first assertion was shown for the case £ — 0). By (12.14j) we have 
A p = 0, so that Supp p C {0} is obvious. We take the Fourier transform of (13.6310 
to find 

{£+ |f1 2 )u-a (f x ■ V s u- u^) + ifip = 0. 

Given 0 G C^M 2 \ {0}) arbitrarily, we set 

poo 

m = / 0(at) e-^ 2 H(O(at) T 0 dt G G 0 °°(M 2 \ {0}), 

Jo 

which solves 

(e + If H 0 + a • V^0 - 0 X ) = 0. 

We thus obtain 

(«, 0) = (% (e + If | 2 ) 0 + a (f X ■ - 0 X )) 

= ((e + If I 2 ) u-a (f x • - n ± ) , 0) 

= ~(i£p,<t>) = 0, 


which completes the proof. 


□ 


The following volume potential (I3.65P is well-defined on account of (I3.23P and 
provides a solution to (I3.7ip for every / G C0°(R 2 ); that is, r o (x, y) is a fundamental 
solution. We also deduce several properties of (I3.65P for later use, including asymp¬ 
totic representation (I3.70p even for less regular /, whose support is not necessarily 
compact but which decays sufficiently fast at infinity. 

Proposition 3.2. Let a G M\ {0}. Suppose 

f G L 1 (M 2 ) fl L°°(R 2 ). (3.64) 

Set 

u(x) = I r a (x,y)f(y)dy, (3.65) 

J R2 

where T a (x,y) is given by (|3.2p . and considerp(x) defined by (13 .5 7p as well. 
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1. The function u(x) is well-defined by (I3.65j) as an element of L'jfifiM. 2 ) (1 S'(R 2 ). 

2. Suppose further that 

/ \x\\f(x)\dx<oo, f(x) = 0(\x\~ 3 (log\x\)~ 1 ) as |x| —» oo. (3.66) 

J R 2 

Then the functions u(x) and p(x) enjoy 

\u(x)\ + |V«(a;)| + |p(x)| = 0(|x| _1 ) as |x| —> oo (3.67) 

with estimate 


sup |x||w(a;)| < C(1 + |a| x ) 

M >p 


(1 + \x\)\f(x)\dx + sup |x| 3 (log |x|)|/(a;)| 

! \x\>p/2 


(3.68) 

for every p > e, where the constant C > 0 is independent of p G [e, oo) and 
a G M \ {0}. Furthermore, we have 

p(x) — / f dy ■—- + 0(|x| -2 ) as hi —>■ oo. (3.69) 

J R2 2n\x\ 2 


3. In addition to (13.641) and (13.661) . assume (12.91) . Then we have 

f 0C^~ 

u(x) — / y L ■ f dy - ——nr + (1 + lal^ 1 ) o(|x| _1 ) as hi — >• oo. (3.70) 
Jr 2 4tt|x| 2 

If in particular the support of f is compact, then the remainder decays like 
o(\x\~ 2 ) m (jsznp . 

4■ Under the conditions (13.641) and (13.661) . the pair {u,p} satisfies 

— A u — a (or 1 • Vu — w 1 ) + Vp = /, cliv u = 0 in M 2 (3.71) 

in the sense of distributions as well as 

(V 2 u, Vp, ■ Vm - G L q (R 2 ) for Vq E (1, oo), (3.72) 

x 2 - ■ Vu G L r (M 2 ) for Vr G (2, oo). (3.73) 

If in addition f G C 00 ^ 2 ), then we have {u,p} G C 0O (R 2 ). 
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Remark 3.2. It is also possible to deduce Wu(x) = 0{\x\~ 2 ) at infinity by use of 
similar estimates of V x T a (x,y), see (|3.3ip . to Proposition \S.1\ (such estimates of 
W x T a (x,y) are not simple consequences of Proposition [701 and one needs further 
several pages). Since slower decay V«(x) = 0(|x| _1 ) in (13.6711 is enough for the 
proof of Theorem, \2.1\ we postpone precise analysis ofW x T a (x,y) until a forthcoming 
paper, in which the external force f = div F with F(x) = 0(|;r|~ 2 ) will be treated by 
using estimates of 'V y T a (x,y). 

Proof of Proposition 1 3 . 21 Let |x| > e, then we take p = 2\x\ in (13.231) to obtain 
[ \T a (x,y)\\f(y)\dy < C(l + |a| _1 )||/|| L oc (R 2 ) |a:| 2 log |x| (|x| > e). 

d\y\<2\x\ 

By (I3.35P we also have 

[ \Ta(x,y)\\f(y)\dy< C f fS + 771732) l/(2/)l d V 

■J\y\>2\x\ J\y\>2\x\ \ 12/1 |®||j/| / 

<C( 1 + |a|^ 1 )||/|| L i (R 2 ) (|x| > e). 

When |x| < e, we similarly use (I3.23p with p = 2e and (I3.35p to find 


|w(x)| < 


+ / < C(\ + |a 


-l 


'|<2e J\y\>2e 


We thus obtain u 6 L“ C (M 2 ) D 2>'(K 2 ). 

We next divide (13.651) into three parts: 

u(x) = Ufix) + U 2 {x) + U 3 (x) 


+ H/||l 1 (k 2 )) (M < e). (3.74) 


+ / + / 

i<M/ 2 J\x\/2<\y\<2\x\ J\y\>2\x\ 

By (I3.33P and (I3.66P we have 

x 1 - f 


T a (x,y)f(y)dy. 


Ufix) = 


47T |cc| 


1<W/ 2 


y ■ f(y) dy + W(x) 


(3.75) 


with 


\W(x)\ < C'|a| _1 |x| _2 
< (27|a|^ 1 |a;|“ 2 


-|<l*l/2 


I f(y)\dy + C\x 


-2 


!) + C\x 


-1 


\y\V(y)\dy 

J\v\<\A / 2 

\y\\f(y)\ dy. 
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The second term of the first line of (13.761) can be estimated even by 


C\x 


-2 


• 1 * 1/2 


(log(e + r)) 1 dr = o(\x\ *) as |x| —> oo. 


Note that this holds true under weaker assumption f(x) = o(|x| 3 ) than (I3.66]) 2 . 
This together with 


x 


47t|x| 2 
implies that 


'\y\> M/ 2 


y L ■ f(y) dy 


< — 


I 3 '! d\y\>\x\/2 


\y\\f(y)\dy = o(\x\ ) 


Ux{x) = 


X 


47t|x| 


y ' f (y) dy + o(\x\ ) as |x| —> oo 


(3.77) 


Let |x| > e, then it follows from (13.2311 with p = 2\x\ and (13.661) that 


\U 2 (x)\ < / \Ta(x,y)\\f(y)\dy 

■' \x\/2<\y\<2\x\ 


< C\x\ 3 (log y) 1 


1 \Ta(x,y)\dy sup |?/| 3 (log \y\)\f(y)\ (3.78) 

\y\<2\x\ \y\>\x\/2 

< C(1 + |a| _1 )|x| _1 sup M 3 (log \y\)\f{y)\ (N>e). 

\y\>\x\/2 

Under stronger assumption (12.91) . we see that U 2 (x) = o(|a;| _1 ) as |x| —* oo. We 
remark that (12.9|) is needed only here. We use (13.351) to find 


\U 3 (x)\<C 


\x\ 1 

+ 


itfi>2|*i \\y\ \ a M 2 


I f(y)\dy 


<C(\x\ 1 + |a| Vl 3 ) / \y\\f{y)\dy = o(\x\ 3 ) 

J\y\>2\x\ 


(3.79) 


as |x| —y oo. We gather (I3.75p . (13.761) . (j3.78j) and (I3.79|) to conclude (j3.68j) for every 
p > e. Then (13.68[) with p — e together with (13.741) for |x| < e yields 


sup(l + |x|)|w(a;)| 

xGR 2 


< C(l + M" 1 ) 


(1 + \x\)\f{x)\dx+ sup (1 + |x| 3 )(log (e + |x|))|/(a;)| 

! xeR 2 


(3.80) 
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Furthermore, we collect (13.771) . ()3.78j) and (13.791) to find the asymptotic representa¬ 
tion (I3.70P as long as (j2.9j) is additionally imposed. If in particular Supp / C B p ( 0) 
for some p > 0, then u(x ) = Ui(x) for |x| > 2 p. In view of the first line of (13.761) . 
we have 

\W(x)\ < C(\a\~ l + p 2 )\x\~ 2 j \f(y)\dy = 0(\x\~ 2 ) as |x| —> oo. 

J\y\<p 

To show the decay of Vu(i), consider 


V(x) := / V x r a {x,y)f{y)dy 


•\<\x\/2 J\ x \/ 2 <\v\< 2 \ x \ J\v\> 2 \ x \ 


= -.V^x) + V 2 {x) + V^x). 


Neglecting the oscillation and using (I3.3ip ~ (|3.32p together with (I3.9p - (l3.10p . we 
deduce 

C\x\~\ \x\>2\y\, ^ 3gl ^ 


|VX,(*,i/)l<| \y\>2\x\. 

Although they are not sharp (Remark 13.2D . they respectively yield 

|hi(x)| < C , |x| _ 1 ||/|| l i ( r2 ) 


and 

/ |y|" 1 |/(y)|dy<C '| a ;|- 2 / |y||/(y)|dy. 

J\y\>2\x\ J\y\>2\x\ 

Let |x| > e and use (13.24p with p = 2\x\ to find 

\V 2 (x)\ < C'|x| -3 (log i^)" 1 f |V x r(x, y)\dy < C|a;| _2 (log -yp) -1 . 

2 J\y\<2\x\ 2 

We thus obtain 

IW(a:)I < ^-r (hi > e). 

fI 

In order to conclude S7u(x) = 0(|;r| _1 ) as |x| —> oo, it suffices to show that 

Wu = V in £>'(M 2 \ B e (0)). (3.82) 

Given <p e (^(K 2 \ B e ( 0)) arbitrarily, we have 


(u, div (p) 


r a (-, y)f{y) dy, div <p) = / (T a (-,y), div <p) f(y)dy, 
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in which the last equality is correct because 


[ [ \T a (x,y)\\f(y)\\dxv (p(x)\dydx <C f oo 

le<\x\<M JR 2 Je<\x\<M FI 


follows from the proof of (13.681) . where Supp ip C Bm( 0) \ B e ( 0). We further obtain 
(u, div ip) = - [ (V x r o (*, y), ip) f (y) dy = -( V ., ip) 


since we have 


\v{x)\ 


/ / |V x r o (z,y)||/(y)||<p(x)| dydx < C / ' r . v 7' dx < oo 

'e<\x\<M JR 2 Je<\x\<M FI 


by computation as above. We are thus led to (13.821) . 
We turn to the decay property of the pressure 


p(x) = 


x 


2n\x\ 


f(y)dy + R(x), 


where the remainder R(x) is divided into three parts: 


R(x) = Ri(x) + R 2 (x) + R3(x) 

1 

2tt 


+ / + / 

<M/ 2 >/|x|/2<|j/|<2|x| J\y\>2\x\ 


x-y 


x 


\x — 2 /| 2 |x| 


f(y ) d v- 


We then observe 


IfliMI < ^ 


3|y| 


\<\x\/2J0 


\x - ty\ 


dt\f(y)\dy< C\x\ - / \y\\f(y)\dy 


and 


H>-( 

7 

J\y-x\<3\x\ 

|xk -i 



1 , 1 
dy + 


V I 1*^1 J\y\<2\x\ 


dy 


as well as 


\Rz(x)\ < 


2vr 


|y|>2|*| 


1 + rh I \f(y)\dy 


\x — y | |x| 


<C\x\~ 2 / \y\\f(y)\dy = o(\x\- 2 ) 

J\y\>2\x\ 
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as |x| —> oo. We thus obtain (13.691) . 

We will show that (13.65|) is a solution to (13.Tip . We use v° and v 1 given by (13.58j) 
and (13.591) . which satisfy (13.601) and (13.611) . respectively, by (13.641) . We set 

poo 

v(x,t) — v°(x,t) + v 1 (x,t), w(x) = / v(x,t)dt. 

Jo 

Since neither u nor w can absolutely converge, we are unable to apply the Fnbini 
theorem directly to them. We will show, nevertheless, that they do converge and 
coincide. Let us employ the centering technique as in (13.26[) . We consider 


u(x) = / Y a (x,y)f(y)dy, 


and 


v(x, t ) = 0(at) J 


p-!/4 1 

K(0(at)x-y,t ) - -^-1 ) f(y ) dy , 


w(x) = / v(x,t)dt, 

Jo 


where T a (x,y) is given by (13.261) . Then both integrals of u and w are absolutely 
convergent over (0, oo) x R 2 with respect to (f, y). In fact, as in (I3.27p . it follows 
from (13.28[) (I3.30|) together with the assumption (13.661) that 


e _1/4t 

K(0(at)x — y-,t) -r——I 


8nt 


1/(2/) I dy dt 


I o Jr 2 


K(x-y,t ) - 


.,—1/44 

-- 1 

8nt 


\f{0(at)y)\dydt 


< C 


(| log |a; - 2/| | + l) 


dy 


j i + \y\ 3 

which is actually convergent. We thus obtain u = w. Since 


POO 


dt 


u -u= 0(at) T e I f(y ) dy = w — w 

Jo Jm . 2 


(3.83) 


and since (j3.83j) does converge by (I3.25p . we eventually conclude that u = w. We 
now show that {u,p} actually satisfies (I3.7ip i in the sense of distributions. Given 
(p G C n °°(R 2 ) arbitrarily, let us consider (u,L- a (p) since we have the adjoint relation 
L_ a = L* a , see (l3ii2l) . Then we hnd 


(u, L_ a ip) = (w, L_ a ip) = / (v(t),L_ a ip) dt, 

Jo 
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in which the Fubini theorem is employed. Note that the argument does not work if 
v is replaced by v. By integration by parts we have 


(u,L_ a p)= / (L a v(t),ip) dt + / (L a (v(t)-v(t)),<p) dt, 


(3.84) 


however, since v — v is independent of x and since 


(L- a ip)(x) dx =—a / ip ± (x)dx , 


we obtain 


(L a (y(t) — v(t)) ,<p) dt — —a(u — u) ± ■ / ip(x) dx 


= a(u — u) ■ / (p ± (x)dx 
J R 2 

= -(u - u, L_ a tp), 

see (13.831) . On the other hand, in view of (13.60[) and (13 .6 1|) and by taking 


(3.85) 


lim {v(t), (p) = 0, lim (v(t) - (/ - Vp), ip) = 0, 

t—> oo t—t 0 


into account, we have 

poo poo 

/ (. L a v(t),<p)dt = - d t {v(t),<p)dt= (f - Vp,(p). 
Jo Jo 

We collect (13.841) . (I3.85P and (I3.86P to obtain 


(3.86) 


(w, L_ a ip) = {f-Vp,tp) 

for all (p G C^°(R 2 ). Since A p = div /, we take the divergence of (I3.71l) 1 to see that 
(div u) G (S'(M 2 ) obeys 

—A(div u) — ax 1 ■ V(div u) — 0 


on account of (12.141) . By Lemma [3.51 div u is a polynomial, however, from (13.671) 
we conclude that div u — 0. Since / G L g (R 2 ) for every q G (l,oo), the result of 
m (see also another proof given by TO implies (13.721) . And then, (13.801) combined 
with (I3.72p especially for leads to (13.731) . Finally, if / G C°°(R 2 ), then 

we put the term x 1 - ■ Vu — u 1 - in the RHS together with such / to use the regularity 
theory of the usual Stokes system. As a consequence, we find {u,p} G C°°(R 2 ). 
This completes the proof. □ 
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For the proof of Theorem 12.21 we also need analysis of the system 


eu — A u — a^x 1 - ■ Vu — u ~ L ) + Vp = /, div u — 0 in K 2 , (3.87) 


where the term eu is introduced in order to control the behavior of solutions at 
infinity. Indeed (13.87j) is the resolvent system, but the only case we are going to 
consider is £ > 0. The velocity part of the associated fundamental solution is given 
by 


T(f\x,y) = / e £t O(at) T K(0(at)x — y,t) dt, (3.88) 

Jo 

while the pressure part is the same, see (13.11) . Of course, (13.881) converges without 
using oscillation, however, what we need is to derive a certain estimate uniformly 
with respect to £ > 0. Therefore, we still use oscillation as well as the centering 
technique. 


Proposition 3.3. Let a e M\ {0}. Suppose f satisfies (13.641) and (13.661) . Set 


u e (x) 


r « ■ ) (x,y)f{y)dy, e > 0, 


(3.89) 


where Fa\x,y) is given by (13.881) . Then u E (x) enjoys (13.681) for every p > e, where 
the constant C > 0 is independent of e > 0 (as well as p G [e, oo) and a G M \ {0 }). 
Furthermore, the pair {u £ ,p} is a solution to (j3.87j) in the sense of distributions, 
where p(x) given by (13.571) . 

Proof. Let m > 0. As in the proof of (I3.7l) - (l3.8p by use of integration by parts, we 
easily find 


< 


'0 


oo Of 

—et+iat p —r 2 /t ^ 
t m 


e 1 e 

C 


+ 




a 2 r 2m 


< 


C 


a r 


2m 


—et+iat 



ds 

^m+1 


dt 


(3.90) 


with some C = C(m) > 0 independent of £ > 0, r > 0 and a£l \ {0}. Owing to 
(I3.90D . we have the similar estimates to (I3.23D . (13.331) and (I3.35j) uniformly in £ > 0; 
namely, there is a constant C > 0 independent of £ > 0 such that 


l<2|x| 


\ 1 ' < a £ \x,y)\dy <C\a\ 1 \x\ 2 + C\x\ 


iT^fa; nil < { ^ + ' |x 

1 a [ ,V)l ~ \ C7|j/|- 1 |ajj + C|aj- 1 | ? /| 


-1UI-2 

j 

-2 


log |x|, |x| > e, 

kl > 212/|, 

\y I > 2\x\. 


(3.91) 
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In fact, it follows from (I3.90p that 



e~ l / At 

8nt 


dt 



with some C > 0 independent of £ > 0, which together with the same computing 
as in the proof of Lemma 13.31 by means of centering technique as in (13.261) yields 
(I3.9ip i . Also, look at the proof of Proposition 13.11 in which oscillation is used in 
(I3.37|) and so on. This time, we employ (I3.90|) to get 



e~ £t O{at) T e~ lxl2/At j 


< 


C 

a\\x 


and so on, where C > 0 is independent of £ > 0. The other estimates without using 
oscillation are obvious. For the purpose here it is enough to split the exponential 
function into two terms rather than (13.361) and '.13. 1 II) since we do not intend to find 
out the leading term. As a consequence, we obtain (I3.91|) 2 . With use of (I3.91|) the 
desired estimate (I3.68P for u e uniformly in £ > 0 is deduced in exactly the same way 
as in the proof of Proposition 13.21 

The proof of the latter assertion is easier than the corresponding part (the 4th 
assertion) of Proposition 13.21 in which we are forced to introduce u. We do not need 
it since u £ itself converges absolutely. Hence, we have 

/•OO 

u £ (x)= / v £ (x,t)dt 

Jo 


with 


v e (x,t) = e £ 0(at) / K{0(at)x — y, t)f(y) dy, 


which satisfies 

d t v £ + (e + L a )v e = 0, v £ (-, 0) = / - Vp 

in M 2 x (0, oo), where L a is given by (13.621) . We thus obtain 


(u £ , (e + L_ a )ip) = (/ - Vp, ip) 

for all ip G Cg^R 2 ). This combined with A p = div / implies div u £ = 0 by Lemma 
13.51 since |Vu e (x)| = 0(|x| _1 ) as |x| —>■ oo, where this decay property is verified 
along the same line as the case £ = 0 by use of (I3.24p and (I3.8ip for V x ri e ^(a;, y) 
without using oscillation. The proof is complete. □ 
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4 


Proof of Theorem 


2.1 


To find the asymptotic representation (12.101) . it would be standard to employ a 
potential representation formula in terms of the fundamental solution T a (x,y) as 
in [9] for the 3D problem, but we have to establish the decay properties (12.7[) in 
advance in order to justify such a formula. This procedure consisting of those 
two steps would be also fine (and actually it works), however, there is another 
way, which is straightforward and leads us directly to (12.101) as well as (12.71) . by 
means of a cut-off technique. We will adopt the latter way to prove Theorem 
12.11 The only disadvantage compared with the former one by use of the potential 
representation formula is that the coefficient of the leading profile needs a bit lengthy 
(but elementary) calculation. 


Proof of Theorem, \2.1\ We use a cut-off technique as mentioned above. In order to 
recover the solenoidal condition by use of the correction term with compact support, 
we first reduce the problem to the one with vanishing flux at the boundary dfl. To 
this end, we fix Xq G int (M 2 \ h2) and introduce the flux carrier 


w{x ) — (3 V ( — log ■ 


-P ~ %o) 


fh r \x — Xq | J 2tt\x — Xq\ 2 ’ 

for given smooth solution {u,p} of (11,5[) . Then we have 


(3 = v ■ u da, 
Jan 


v ■ w da = (3, 


'an 


and 


div w = 0, Aw = 0, (x — xq) ± ■ Viv = w J 


V J w(x ) = ) + 0(|x|- (2+j) ) 


27t|x| 


in R 2 \ {iCo} 

0 = o, i) 


as |x| —» oo. So the pair 


fulfills (11.5)) subject to 


u = u — w, p = p — ax 0 ■ w 


v ■ u da = 0, 


i an 


where we note the relation 


4W' «'i = (J: log |. r 1 — J 


(4.1) 

(4.2) 


(4.3) 


= Xq -Vw k (k — 1,2). (4.4) 
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We fix R > 1 such that R 2 \ C B R { 0). Let i/j E C™(B 3R ( 0); [0,1]) be a cut-off 
function satisfying ip{x) = 1 for |x| < 2 R. By using the Bogovskii operator B in the 
annulus 

A = {x E R 2 ; R < |x| < 3-R}, 

see a, 0 and [T9] , we set 

v — (1 — ifj)u + B[u ■ V-0], q = (1 — i>)p- 

It should be noted that f A u ■ 'Vijjdx = 0 follows from (14.3k Then the pair {v,q} 
obeys 

— Av — a (or 1 ■ Vu — v 1 ) + Vq = g + (1 — div v — 0 in R 2 (4.5) 

for some function g E C£°( R 2 ) whose support is a compact set in A. Here, we do 
not need any explicit form of g ; in fact, the important quantity (14.81) below can 
be calculated only by taking account of the structure of the equation (14.51) . that 
is, div S(v,q ) = —g — (1 — ^)/, see (12.41) . When u(x) = o(|x|), it is obvious that 
v E iS^R 2 ). Under the alternative assumption Vu E L r (Q\B R (0)) for some r < oo, 
we have Vu G 5'(R 2 ), which implies v E iS'(R 2 ) by [7, Proposition 1.2.1], Going 
back to (14.5p . we observe Vq E 5'(R 2 ) and thereby q E 5'(R 2 ), too. Proposition 13.21 
together with Lemma 13.51 concludes that 



(4.6) 


with some polynomials P v and P q , however, it turns out from (13.671) and from either 
Vu G L r (R 2 ) with some r E (l,oo) or v(x) = o(|x|) that P v must be a constant 
vector Uoq. Thus we have 



u(x) = w(x)+ r a (x,y){g+{l-'ip)f}(y)dy + u 00 (|x| > 3R), (4.7) 


from which combined with (13.701) and (14.21) we obtain (12.101) under the additional 
condition (12.91) as well as (12.71) -^ where the coefficient 


a= f y^ ■ {g + (1 -ip)f}(y) dy = - [ y ± ■ div S(v, q) dy (4.8) 



is computed as follows. 
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Set 


«(p) : = - / y ± ■ div S(v,q)dy (p > 3/2). 
J\y\<p 


In view of (12.41) we have the relation 


to find 


div (y ± ■ S(v,q )) = ^d k [(y^jS^v, g)] 

j,k 

= y 1 - • div S(v, q) - S 12 (v , q) + S 2 i(v, q) 
= y L ■ div S(v,q) — 2 ay ■ v 


a (p) = ~ V ± ■ ( S(u,p) - ) da - 2a / y-vdy. 
J\y\=p \ P/ J\y\<P 


(4.9) 


Since div S(u,p) = —/ in fl, it follows from (j4.9j) in which v is replaced by u that 

a{p)= j y ± ■ (S(u,p)u) da+ 2a y-(u-v)dy+ y L ■ f dy. 

J dn Jfip Jfip 

We are going to compute 


y ± ■ (S(u,p)v) da 


I dn 


>dn 


We will show that 


= / y L ■ {(T(u,p) + au ® y L ) u] da 
Jdn 

— / y^ ■ ((Dw)v) da + a / (y L ■ v) (x,} • w) da 

Jdn Jdn 

— a y L ■ {(w ® y^v} da — a y ± ■ {( y ~ L ® u)v} da 

Jdn Jon 

y ± ■ { (' T(u,p ) + au® y 1 -) v) da + J\ + J 2 + -h + J±- 


J] — 0, J 2 4” J3 — 0, J 4 + 2a / y ■ (u — v) dy = 0, 

J fin 


which concludes 


<*(p)= y ± ■ {(T(u,p) + au®y ± ) v] da + y ± -fdy. 
Jdn Jn B 


36 




Letting p —> oo leads us to 


a 


= / • {(T(w,p)+ 0 ^®^) z/} da+ / y ± -fdy. 

Jan Jn 

In fact, we observe 

2 a y-(u — v)dy — a {^u — B[u ■ V-0]} • V|j/| 2 dy 

Jn p Jn p 

— a f div [ \y\ 2 {ipu — B[u ■ Vgb]} ] dy 


(4.10) 


and 


= a \y\ 2 (y ■ u) dcr = — J 4 
Jan 


J2 + J3 = —a / (y -1 • is) (y — ^o)" 1 • w da = 0 
Jan 


on account of (y — ^o)" 1 • w(y) = 0. We take account of (Vtc) 2 = Vw and A w = 0 
in M 2 \ {xq} to see that 


J\ = — 2 / (y — xq) l ■ (v ■ X7w) da — 2 xq ■ / v-Wwda 
Jan Jan 

— —2 f (y — xo)" 1 • (z' • Vro) da + 2xq • f ——— 

*/,9Q J\y—x o|=e ^ 


■ Vw dcr 


where £ > 0 is taken in such a way that B e (x 0 ) C int (IR 2 \ Q). Using the explicit 
representation 


Vw(y) = 


P ( I 


2 tt \\y-x 0 


we find 


'an 


[ 2(y - x 0 ) <g> (y 

Xq\ 2 

|l/-a;o| 4 

~P 

f (jj - xo) 1 - ■ v 

2vr 

Jan \y~x o| 2 

J_ 

/ div 

2 7T „ 

/R 2 \(nuB e (a:o)) 


da 


■dy = 0 


and 


y~x 0 P 

-- vw da = 


'\y-xo\=e e 


27T£ 3 


’\y-x o|=e 


(y - ^ 0 ) do- = 0 
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which implies that J\ = 0. We thus obtain (14.101) . 

Concerning the pressure, it follows from (j4.4j) and (14.61) that 

Vp = aXo • Vw + V / Q(x-y)-{g+(l-i())f}(y)dy + VP q (|x| > 3R). 

Jr 2 

By (13.721) together with (14.21) we know V(p — P q ) G 17(0) for every r G (l,oo). 
Since Aw = 0 in M 2 \ {xo}, we obtain from (14.71) 

Au = aI T a (x,y){g+ (1 - i>)f}(y) dy (|x| > 3R), 

Jr 2 

so that A u G 1/(0) for every r G (1, oo) on account of (13.72p . In addition, we also 
have 

x x • Vu = x x • \7w + x L ■ V / r a (x, y){g + (1 — dy (|x| > 3R). 

Jr 2 


It thus follows from (13.731) and (14.21) that x 1 - ■ V« G L r (h2) for every r G (2,oo). 
Taking those as well as (EZDi into account, we go back to (1 1.5 1) and let |x| —> oo to 
find that VP, = — au^. This implies that 

p = clxq ■ w + f Q( X ~y)-{g + {l-il>)f}{y)dy-aui'X + p eo (|x| > 3 R) 
Jr 2 


for some constant p^. By (I3.67[) together with (14. 2 p we obtain (I2.7L . We also use 
(I3.69P and carry out a bit computation to obtain 


p(x) + au^ ■ x - poo = / [T(u,p) + a [u ® y L - ® u)} vdcr y 


'an 


f dy - /3aXo 


x 


(4.11) 


27T X 


+ 0( |x|' 2 ) 


as |x| —> oo. We stop further computation of the coefficient, however, we will recall 
(14. lip in Theorem 12.21 in which the coefficient is much simplified. 

Once we have fine decay properties (12.71) . we are able to justify the energy relation 
(12.81) . We first verify (12.5p for smooth vector fields u, v G H^JQ) without assuming 
their decay properties at infinity. For each p > 0 large enough we have 



[(x -1 • V« — u L ) • v + u ■ (x -1 ■ Vu — u -1 )] dx 


div [x ± (m ■ v)\dx 


{u ■ x^){u ■ v) da 


'an 


(4.12) 
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since f , = 0. Letting p —> oo, we obtain (12.5[) . Now, given smooth solution 

{w,p} G x L 2 oc (h2), we use the constants {«oo,Poo} found above and set 

m*(x) = u(x) - Moo, p*{x) = p(x) + au ^ • x - Poo, 


which satisfy 


—Am* — a (or 1 • Vm* — u^) + Vp* = /, div m* = 0 in 0. 


We multiply u*, perform integration by parts over Q p and use (I4.12|) to find the 
following two equalities, in which the relation div T(u*,p*) = Am* — Vp* is used for 
the latter: 


1 

2 



| Vm*| 2 dx 



| Dm* | 2 dx 



pj)v) ■ u*da + /, 



(T(m*,p*)z/) 


• M* d(T + /, 


where the common term / is given by 


/ 


• 2T~)|m *| 2 dcr + / f-u*dx. 


i an 


Note that both T(u*,p*)z/ and (Vm* — p*I)z/ are understood as the normal trace 
being in H~ 1 ^ 2 (dD, p ). In view of (14.71) . we see from (13.671) and (14.21) that 


Vm*(i) = Vm(i) = 0(|x| x ) as |x| —>■ oo. 


This together with (12.71) implies that 


lim / 

^°° d|x|=p 



{y ■ M*)p* ) da 


0, 


and that 


lim / (T(m*,p*)^) • M*d(T = 0 

^ >:x; J x| : 


On the other hand, we know that / • m* G L 1 (O) by (12.6jl and (12.71) together with 
m G C LlM for all s < oo. We thus obtain not only Vm G L 2 (Q) but 

(12.81) . This completes the proof. □ 
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5 


Proof of Theorem 


2.2 


Proof of Theorem, \2.S\ We begin with the proof of uniqueness. Suppose {u,p} is a 
solution in the sense of distributions to ca with / = 0 subject to u = 0 on <912 
and {u,p} —> {0,0} as |x| —> oo within the class Vw G L 2 (f2), {u,p} G LlM- We 
put the term x^ ■ Vu — u ^ in the RHS and use the interior regularity theory of the 
usual Stokes system to show that u G H^ 1 (12), p G Hf oc (il) for every integer k > 1 
by bootstrapping argument; hence, m,})G C' 0O (n). By Theorem 12.11 we have (12.81) . 
in which the RHS vanishes. So, u is the rigid motion, but u — 0 on account of the 
boundary condition. Going back to (j 1.5 [) (with / = 0), we have Vp = 0, which 
together with p —> 0 at infinity yields p = 0. This proves the uniqueness. 

We turn to the existence. It is easy to find a solution with Vu G L 2 (12) by 
following the method of Leray, but one cannot exclude a constant vector u 0 0 at 
infinity even if applying Theorem 12.11 To get around this difficulty, we will adopt an 
approximation procedure specified below which brings regularizing effect at infinity. 
We take the auxiliary function 

C'(M) + C(M)j 

where ( G C°°([0, oo); [0,1]) satisfies ((r) — l (r < R) and C(r) = 0 (r > 2 R), where 
R > 1 is fixed such that M 2 \ H C Rr( 0). Then we have 

w\on = ax' 1 , div w = 0, x 1 ■ \7w — w ± = div (w <g) x 1 - — x 1 - ® w) = 0. 


(ax -1 ) (5.1) 



We will find a solution of the form u = u + w, where u should obey 

{ —A u — a (x -1 ■ Vu — u 2 -) + X7p = f + Aw, div u = 0 in H, 
u\an = 0, u —y 0 as |x| —> oo. 


(5.2) 


As Finn and Smith performed in their paper pjjj on the Oseen system (see also 
Galdi [19;, Section VII.5]), for e G (0,1), let us consider the approximate problem 

{ eu £ — A u e — a (x -1 ■ Vu e — u^r) + S7p e = f + Aw, div u £ = 0 in Q, 
u £ \dn = 0, u £ —y 0 as |x| —> oo. 

(5.3) 

By (^^.(H) we denote the class of all solenoidal vector fields being in (^(fl). Let 
be the completion of 0^(12) in // 1 (f2). In a usual way (see, for instance, 
the proof of Lemma 5.3 of [27] . in which the problem in each bounded domain Q p 
is first solved by means of the Lax-Milgram theorem and then the limit p —» oo is 
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considered by using a priori estimate uniformly in p), one can find u £ G H ( \ ^(Q) 
which satisfies 

e llwlll 2 ^) + 2 Il^ 7 ' u ell2 2 (0) - 2 + ^Hli 2 (n) (5-4) 

and 

e(« e , ¥>) + (Vu e , V(/p) - a {x^ ■ Vu £ - uf, p) = (/ + Aw, 99 ) 

for all y? G Co^-(^)- We choose p £ G L z 2 oc (h2) such that J n p £ dx = 0 and that the 
pair {u £ ,p £ } satisfies (15.31) { in the weak sense. Since / + Aw G C°°(f2), the interior 
regularity theory of the usual Stokes system implies that u £j p £ G C°°(r2). 

As in the proof of Theorem 12.11 we take the same cut-off function 0 together 
with the Bogovskii operator B in the annulus A = {x G M 2 ; R < |x| < 3R} and set 

v £ = (1 - ip)u e + B[u £ ■ V0], q e = (1 - ip)p e . 

Then the pair {u e , q £ } obeys 


ev £ - Av £ - a (x x • Vn £ - v^) + Vq £ = g E + (1 - i/j)f in M 2 , 

div v £ = 0 in M 2 , 


where 


g £ = eB[u £ • V0] + 2V0 ■ Vu £ + (A-?/; + aar 1 • Vi/))u e — A B[u e ■ V0] 

— ax -1 • VB[u e ■ V-0] + aB[u £ ■ V^] ± — (Vip)p e . 

Here, note that (1 — -0)Aw = 0 since -0 = 1 (|x| < 2R) and since Aw = 0 (|x| > 2 R). 
We use the fundamental solution (13.881) for the system (13.87(1 . Then, by Proposition 
13.31 with p = 6 R and Lemma 13.51 we find 


v e {x) 


r i e) 0A y) {9s + (1 - 0)/} ( y ) dy 


subject to 


sup |x||n e (a;)| < C (1 + |a| x ) 

|x|>6ii 


[ ( 1 + kl)|{^ + (l-V')/}(^)l dx 
Jr 2 


+ sup |x| 3 (log |x|)|/(a:)| 

|x|>3i? 


(5.5) 


with C > 0 independent of e. Here, the point is that a constant vector can be 
excluded since u £ G L 2 (H). 
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By f n , iR Pe dx = 0 and (15.31) , we have 

\\Pe\\L 2 (n 3 R) ^ C'fl||Vp.|U-i (n3ll ) < C R (||M £ ||^i(f 2 3H ) + \\F + Vw|| L 2 (n3i?) ) , 

where H^ 1 (Q 3 r) := This together with (15.41) and the estimate of the 

Bogovskii operator (Bl. ® and [19j ) lead us to 

[ \ 9 e(,y)\dy < 2 V 2 R\\g £ \\ L 2 {A) 

J A 

< C R (||n £ ||//i(o 3fl ) + ||p £ ||L 2 (n 3i? )) 

< Cr (11 Vrte11 r 2 (q 3R ) + \\F\\L 2 (n 3R ) + M) 

< C7fl(||F|| La(n ) + |a|), 

which combined with (15.5|) implies that 

K(z)| = |u e (a:)| < C( 1 + |a| _ 1 )(|a| + ||F|| L 2 (n) + [/]) |x | _1 (|x| > 6 R), (5.6) 

where 

[/] := j \x\\f(x)\dx+ sup |x| 3 (log|x|)|/(x)| 

J\x\>2R hl>3K 

and C = C(R ) > 0 is independent of e E (0,1). By (j5.4j) we have 

llwll l 2 .°°(o 6J? ) < C\\u £ \\ L 2 (p. 6R ) < C R \\Vu e \\ L 2 {n6R) < C R (\\F\\ L 2(ty + |a|), 
which together with (15.6|) yields 

u £ E L 2 '°°(Q), ||m £ ||l 2 .°o(o) < C {l + M + (1 + M *) (II-^IIl 2 ^) + [/])} 

with C = C(R ) > 0 independent of £ E (0,1). Hence, there is u E L 2 ’°°(Q) (1 Hl oc (Q) 
with Vw E L 2 (Q) such that, as e —* 0 along a subsequence, 

u £ —» u weakly-star in L 2 ’°°(Q), Wu £ —> Vn weakly in L 2 (fi), 

u £ —y u weakly in if 1 (0 /0 ), u £ —>■ u strongly in L 2 (Q p ), 

for every p > R and, thereby, 

(Vn, V<p) — a (x 3- • Vm — hr 1 , <p) — (/ + Aw, p>) 

holds for all <p E as well as div u = 0. We fix p and use the trace inequality 

K — u \\L 2 (dvip) A C\\u e — u\\ £ 2 (Q p )ll' u e — u \\m(p, p ) 
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to see that u\gn = 0. Since A u + a (or 1 • Wu — u 1 -) + / + A w G H 1 (Q p ) for 
every p > R, we find an associated pressure p G L 2 OC (0) such that the pair {u, p} 
solves (15.21) , in the weak sense. The interior regularity theory of the Stokes system 
concludes that {u, p} is smooth and, therefore, so is u := u + w. Both estimates in 
(12.15P are obvious. 

Let us apply Theorem 12.11 to {u,p} with Vn G L 2 (0) as well as u G H^ oc (Q), 
p G L 2 oc (Q). Since u G L 2 , 00 (f2), we have all the properties in this theorem with 
Uoo = 0. We denote p — p 0 0 by the same symbol p so that {u,p} is the desired 
solution. By w|an = ax 1 - we have 

P= v-udcr = 0, / y L ■ {(w ® da y = a / {y ■ y^\y\ 2 da y = 0, 

Jon Jon Jan 

and thereby (12.101) implies (12.161) . Finally, asymptotic representation of the pressure 
is given by (14.111) . in which {woo,Poo} = {0,0}. Since j3 — 0 and u\gn = ay 1 -, we 
conclude (12.171) . The proof is complete. □ 
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